Can a local repulsive potential trap an electron? 
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We study the classical dynamics of a charged particle in 
two dimensions, under the influence of a perpendicular mag- 
netic and an in-plane electric fleld. We prove the surprising 
fact that there is a finite region in phase space that corre- 
sponds to the otherwise drifting particle being trapped by a 
local repulsive potential. Our result is a direct consequence 
of KAM-theory and, in particular, of Moser's theorem. We 
illustrate it by numerical phase portraits and by an analytic 
approximation to invariant curves. 



The detailed dynamics of electrons in two dimensions 
(2D), in electromagnetic fields and with localized scatter- 
ers, is important for several aspects of the quantum Hall 
effect. Classical dynamics is more relevant in this context 
than one might think, see the recent review by Trugman 
. In addition, classical kinetic theory for such systems 
is surprisingly subtle ||^, and requires for its foundation 
an understanding of the underlying dynamics. In this 
letter we focus on the classical dynamics of an "electron" 
(i.e., a particle with charge — e) confined to the xy— plane. 
Perpendicular to these two dimensions there is a constant 
magnetic field B — Be^. In addition, there is a constant 
in-plane electric field £ = £ey. When no forces act on 
the electron except those stemming from these fields, the 
nature of the motion is well known. With £ ~ 0, B ^ 
the electron with energy E — ^mv^ moves in a circu- 
lar orbit with cyclotron radius R = v/ojc and cyclotron 
frequency lOc = eBjm. Addition of a finite electric field 
results in the center of this orbit acquiring a constant 
drift velocity Vd = {£ 'x. B) / B'^ = {£/B)ex- The question 
we address in this letter is the following: With both £- 
and ;B- field present, is there a finite region in phase space 
which corresponds to the otherwise drifting electron be- 
ing trapped by a single scatterer, i.e., by a local repulsive 
potential? (For concreteness, we let the repulsive poten- 
tial be a hard disk of radius a.) This question arises, for 
example, when one tries to construct a kinetic theory for 
the Lorentz gas in 2D, with both electromagnetic fields 
present. 

When £■ = 0, B 7^ all electrons are, in a sense. 



trapped by the magnetic field, and the hard disk simply 
modifies the cyclotron orbits into orbits skipping around 
the periphery of the disk (see for instance 1^). On the 
other hand, when £ ^ B = 0, one or more collisions 
with the disk will only temporarily perturb the dynam- 
ics determined by a constant acceleration in the negative 
electric field direction -Gy. With both fields present, the 
question is whether an electron initially colliding with 
the disk, will ultimately miss it, due to the drift caused 
by the electric field. 

There are two dimensionless parameters in the prob- 
lem, the ratio of the cyclotron radius and that of the disk, 
r = R/a, and the ratio of the displacement, Vd ■ 2n/LUc, 
during one cyclotron period and the radius of the disk, 
e = (27rm/ae) • {£/B^). Clearly, the answer to our ques- 
tion whether the disk can trap the electron when drift is 
included, must depend on where in this two-dimensional 
parameter space it is posed. Our aim in this letter is not 
to give a complete answer to the question, with parameter 
regions and corresponding measures precisely delineated. 
Our more modest aim is to demonstrate the surprising 
qualitative fact that a repulsive potential can, indeed, 
trap a charged particle in crossed electromagnetic fields 
in 2D. 

We first emphasize that an electron drifting towards 
the hard disk from "infinity" cannot be trapped. The 
main elements of the proof are: To the first collision with 
the disk is associated a finite measure in phase space. 
Consider the subset of those incoming trajectories (if they 
exist) that will hit the disk at least N times. Since classi- 
cal dynamics is measure preserving, the regions in phase 
space associated with each of the N collisions have the 
same measure. Reversibility (invariance under reversal of 
time and magnetic field) implies that they do not overlap. 
Since the total phase space associated with an electron 
colliding with the disk is finite, it follows that the mea- 
sure of the trajectories coming from infinity and hitting 
the scatterer at least N times has to converge to zero 
as TV ^ oo. (It may, of course, vanish for a finite TV.) 
Trapping is, therefore, only possible for electrons that are 
initially close to the scatterer. 

Naively one might think that conditions are more fa- 
vorable for trapping when r ^ 1. If this were so, not 
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much room would be left for trapping: When r ^ 1, 
the perimeter of the disk can, on the scale set by the cy- 
clotron radius, be considered a straight line. An electron 
skipping upwards along the right side of the disk will, 
unless the curvature of the disk makes itself felt in time, 
sooner or later miss the disk and drift away from it. 

Surprisingly, it is the opposite regime, when r ^ 1, 
which is more favorable for trapping. Here we can, on 
the scale of the radius of the disk, consider the parts of 
the electron orbit close to the disk as straight lines. Two 
subsequent collisions are shown in Fig.|l|. The direction 
of the motion just prior to collision no. n is defined by 
the angle 4>n- To the extent that the orbit locally can 
be considered to consist of straight lines, the scattering 
angle of collision no. n is given by ipn — 4'n+i — 4>n- Due 
to the electric field, the incoming line heading for collision 
no. (n + 1) is shifted by the amount ea in the positive x- 
direction with respect to the outgoing line from collision 
no. n. In terms of the angles (/)„ and the dimensionless 
impact parameters /?„ = 6„/a, inspection of Fig.|l| gives 
the map, 

(l>n+i = 0„ + TT - 2 sin"^ /3„ 

Pn+l ^ Pn-e sin (j)n+l- (1) 

In contrast to standard scattering conventions in three 
dimensions, inherent in the present map is the following 
choice: The impact parameter (3 lives on the interval 
(— 1, 1) and the scattering angle -0 on the interval [0, 27r), 
with /3 — > — 1 corresponding to "0 — > 27r. 

One can show that the map resulting from the full 
dynamics for arbitrary B reduces to the map (|^) when 
the time interval between successive collisions is approx- 
imated by the cyclotron period 2it /ujc- This implies that 
corrections to are uniformly bounded by a constant 
times a/R = r^^ . Consequently, the map (P is reliable 
when r ^ 1. 

The map has several interesting properties. First of 
all, since its Jacobian equals unity, it is area-preserving. 
Thus, our weak S-field approximation, basic to (^), has 
exactly preserved this general property of Hamiltonian 
maps. Next, for small e, (Q) is a small perturbation of 
the integrable map 

(i>n+l = 4>n + ^{Pn) 

dn+l^Pn, (2) 

where fi(/3) = tt — 2 sin^^ j3 describes scattering off a hard 
disk. From (H) we recover the well-known fact that the 
scattering angle = 4'n+i — 0n remains the same in 
subsequent collisions, when the electric field vanishes 
Furthermore, ^{(3) is a monotonic function of its argu- 
ment and, as a result, our map (Q) belongs to the category 
of twist maps [Q , for which a number of properties have 
been rigorously established. 

However, the question of principal interest from our 
point of view is the existence of trapped orbits, and the 



phase space associated with them, when e > 0. The 
function ri(/3) is undefined when \(3\ > 1. Physically, this 
corresponds to the electron missing the disk. The regions 
of first and last collision are represented in phase space by 
areas of order e near j3 = +1 and /3 = — 1. The nontrivial 
question is, therefore, whether there is a finite region of 
phase space associated with orbits that never enter the 
escape region. 

For small enough e, Kolmogorov-Arnold-Moser 
(KAM) theory enables us to give a positive answer to this 
question. Specifically, in our case of a two-dimensional 
map, we will use Moser's theorem |^ which proves the ex- 
istence of invariant curves in phase space for small pertur- 
bations, and for some initial conditions. The existence of 
two distinct invariant curves winding around phase space, 
together with area preservation, shows the finite region 
enclosed between these two curves to be an invariant set, 
corresponding to trapped orbits. It is straightforward to 
verify that the map (or any small perturbation of it 
that remains area preserving) satisfies the hypotheses of 
Moser's theorem, provided that we restrict phase space 
to a smaller strip < c < 1, in order to avoid the escape 
regions and the singularities of 17. The important point is 
that this theorem does not require the entire strip to be 
invariant under the map. We have thus proved the exis- 
tence, for small electric and magnetic fields, of invariant 
curves in phase space, and hence of trapped orbits with 
positive measure. 

This rigorous result is illustrated numerically by fig- 
ures 1^ and ^ where we show two different phase por- 
traits of the map. When e — 0.2, many invariant curves 
are present. We can divide them into two classes: those 
winding around phase space, which we study in this let- 
ter, and those forming "elliptic islands" around periodic 
orbits. When e = 0.4, the structure of phase space has 
become more complex. Invariant curves winding around 
phase space are now scarce and, numerically, all of them 
seem to be destroyed between e — 0.4 and e = 0.5. On 
the other hand, elliptic islands, which are numerous in 
Fig.^, survive at much stronger electric fields (certainly 
up to e 1.1). However, most of the chaotic trajectories, 
which remained bounded for small e, diffuse until they 
reach the escape region. 

As an analytic illustration of our rigorous result, we 
have also constructed a perturbation scheme around (H), 
close in spirit to that used in the proof of Moser's theorem 
We start from the Ansatz 

4>n = C,+nijj + eui{n; (,lu) + e^U2{n] C,,bj) + ■ ■ ■ (3) 
V'n = 0„+i - 0n = + ewi (ti; C, lu) + £^W2 (n; C, w) H , 

with Viiji) = Ui{n+1) — Ui{n). Both the "effective initial 
condition" C, and the "angular frequency" oj depend on e 
and on the true initial conditions (0Oj ''Po)- Setting n — Q 
in (^) gives the inverse relations = 4>[){C,,oj;£), ipo = 
V'o(Ci'^;£) via the functions Ui{0-X,uj), Ui(0;C,w). We 
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determine Ui{n) and Vi{n) by inserting the Ansatz (|^) 
into the map ([^) and doing straightforward perturbation 
theory, treating C and lo as constants in the process. We 
fix the "initial values" Wi(0) and Vi{0) by requiring that 
the oscillating functions Ui{n) and Vi(n) average to zero 
for large n. (This is indeed necessary for the Ansatz 
to be consistent.) Convergence of the scheme is assured 
if a;/27r is a Diophantine number (loosely speaking: 
an irrational number which is hard to approximate by a 
sequence of rationals). 

To zeroth order the result is simple, C = 0Oi = "001 
reproducing the skipping orbits in zero electric field. To 
first order, perturbation theory gives 



vi{n) 



sm\uj/2) 



{cos[C + w/2] 



.s[C+(n+iV]}+wi(0). 



(4) 



In (^) we determine vi (0) by requiring vi (n) to average to 
zero for large n. We can then calculate ui{n) by direct 
summation, again fixing ui(0) by the stipulation that 
ui{n) averages to zero. The results are, 



, cos(C + w/2) sinC 

vi(0) — ^ : UilO) = 5 . 

^ ' sin2(tj/2) ^ ' 2sin3(w/2) 



(5) 



In this way one can continue. We are primarily inter- 
ested in the frequency w(e). To second order, we deter- 
mine W2(0) in the same manner, and contruct the inverse 
relation from to 0{e^) as 



cos(C + w/2) 2 
^ sin2(w/2) [ 



C0S(2C -I- Ld) 



sin'^(tj/2)cos(u;/2) 



(6) 



Inversion of this relation (with appeal to the first order 
result for C,) to 0(2^) finally yields 



cos 



io + ^-0/2) 



sin2(Vjo/2) 
cos(?/'o/2)[cos(20o 



— e 



■ V'o) 



cos(2(/)o + iPq) 
8sin^('0o/2)cos(-0o/2) 
f3] 



4sin^(V'o/2) 



(7) 



Note that in the second order term of (^) , there appears a 
"small denominator", cos(ti;/2), which vanishes at uj = n, 
i.e., (3 = (where an orbit of period two exists). Higher 
order terms will generate a proliferation of other small 
denominators. This denominator appearing in second 
order is the first symptom indicating that convergence 
of the scheme requires lu/2'k to be a Diophantine num- 
ber. To see that the solution indeed describes an invari- 
ant curve, it is sufficient to note that, by construction, 
Vi{n; Lu) = Vi{0; ( + ncu, lu), with a similar equation for 
Ui. Thus, for given e and lu, (jjn and ipn depend only on 
the combination Q + mo. Elimination of this parameter 
and use of (^ yields an explicit equation for the invari- 
ant curve, fin = ^~^{'>Pn) = /(0n;w,e)- As stated above, 
once two distinct invariant curves have been constructed, 



they delimit an invariant region of positive measure cor- 
responding to trapped electrons. 

For Diophantine when the perturbation scheme 

implied by (^ converges, by taking the average one con- 
cludes that (V'n) = w. In Fig.^ we show a comparison 
of the scattering angle averaged over 1 000 000 collisions, 
and our analytic perturbation result for the frequency, 
with the initial conditions (0o = 0.3, ■00 = 2.0). For 
small e the results are in very good agreement, despite 
the fact that w takes both rational and irrational values. 
For larger e, one can observe the breakdown of perturba- 
tion theory, due to the phenomenon of "mode locking": 
On the interval of roughly 0.20 < e < 0.45 the electron 
is trapped in the vicinity of a periodic orbit with period 
3 and, as a consequence, (V'n)/27r = 1/3, a rational num- 
ber. 

In conclusion, we have shown that for sufficiently small 
electric and magnetic fields, bound states associated with 
a hard disk scatterer constitute a set of positive measure. 
Indeed, both fields add small and smooth perturbations 
to the integrable map (||), so that KAM theory applies. 
By investigating the stability of periodic orbits, it is in 
fact possible to show that trapping occurs also for mod- 
erate values of the magnetic field. Note, finally, that 
generalization from a hard disk to an arbitrary, rotation- 
ally invariant repulsive potential of strictly finite range, 
leaves the scattering function 51 (/3) monotonic. In other 
words, the corresponding map is still a twist map, and 
all qualitative conclusions remain unchanged. In short, 
we have accomplished our goal and, within classical me- 
chanics, answered the question: Yes, a local repulsive po- 
tential can trap an electron! With quantum mechanics, 
the question remains open. 

Acknowledgements. We would like to thank Professors 
J. S. H0ye and H. Kunz for important remarks. N.B., 
E.H.H. and J. P. are grateful to Professor J. -P. Hansen 
and the Ecole Normale Superieure de Lyon, for their hos- 
pitality. Finally, N.B. ackowledges financial support from 
the Swiss National Fund, J. P., A.H. and E.H.H. from 
The Norwegian Research Council and J. P. from KBN 
(The Polish Committee for Scientific Research), grant 
no. 2 P302 074 07. 



[1] S. A. Trugman, Physica D 83, 271 (1995). 

[2] A. V. Bobylev, Frank A. Maa0, Alex Hansen and E. H. 
Hauge, Phys.Rev.Lett. 75, 197 (1995); J.Stat.Phys. (in 
press); J. Piasecki, Alex Hansen and E. H. Hauge, to be 
published. 

[3] N. Berglund and H. Kunz, J. Stat. Phys. 83, 81 (1996). 
[4] J. D. Meiss, Rev. Mod. Phys. 64, 795 (1992). 
[5] J. Moser, Nachr. Akad. Wiss., Gottingen, Math. Phys. 
KL, 1 (1962); Stable and Random Motions in Dynamical 



3 



Systems (Princeton University Press, New Jersey, 1973). 
[6] Some readers may find the following motivation for the 
Ansatz useful: An invariant curve has a parametric equa- 
tion of the form = ^ + £U{^,u),e), ip = to + £V{^,u),e). 
The parameter ^ grows proportionally to time with "an- 
gular frequency" ui: ^„ = ^ + nui. 

FIG. 1. Geometry of two successive collisions with the 
hard disk. 

FIG. 2. Phase portrait of the mapping for e — 0.2. The 
escape regions are recognized as those close to the upper and 
lower part of the figure. Note that orbits of period 2, 3 and 
4 are easily located. Between the periodic orbits are bands of 
invariant curves, described by the analytic theory of the text. 

FIG. 3. Same as Fig.2, but for e = 0.4. The escape regions 

have grown in size, more periodic orbits can be recognized, 
and the regions with invariant curves have shrunk consider- 
ably. In addition, chaotic regions (which remain trapped!) 
have become manifest. 

FIG. 4. Comparison between the "frequency" a;/27r as 
computed from our analytic approximations to first and sec- 
ond order, and (tpn)/2Tr obtained from the simulation, where 
the average has been taken over 1 000 000 collisions, starting 
from the initial condition (po — 0.3, i/'o = 2.0. Note, with ref- 
erence to Fig.2, that the electron is mode locked to a period 
3 orbit, {ip„)/2-!T = 1/3, from e w 0.2 to e w 0.45. 
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